Spin excitations in a monolayer scanned by a magnetic tip 
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Abstract. - Energy dissipation via spin excitations is investigated for a hard ferromagnetic tip 
scanning a soft magnetic monolayer. We use the classical Heisenberg model with Landau-Lifshitz- 
Gilbert (LLG) dynamics including a stochastic field representing finite temperatures. The friction 
force depends linearly on the velocity (provided it is small enough) for all temperatures. For 
low temperatures, the corresponding friction coefficient is proportional to the phenomenological 
damping constant of the LLG equation. This dependence is lost at high temperatures, where the 
friction coefficient decreases exponentially. These findings can be explained by properties of the 
spin polarisation cloud dragged along with the tip. 
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, Introduction. — While on the macroscopic scale the 
phenomenology of friction is well known, several new as- 
pects are currently being investigated on the micron and 
nanometer scale [1,2]. During the last two decades, the 
r,esearch on microscopic friction phenomena has advanced 
enormously, thanks to the development of Atomic Force 
Microscopy (AFM, [3]), which allows to measure energy 
dissipation caused by relative motion of a tip with respect 
tio a substrate. 

' Recently the contribution of magnetic degrees of free- 
dom to energy dissipation processes has attracted increas- 
ing interest [4-8]. Today, magnetic materials can be con- 
trolled down to the nanometer scale. New developments in 
the data storage industry, spintronics and quantum com- 
puting require a better understanding of tribological phe- 
nomena in magnetic systems. For example, the reduction 
of heat generation in reading heads of hard disks which 
work at nanometer distances is an important issue, as heat 
can cause data loss. 

Magnetic Force Microscopy (MFM), where both tip and 
surface are magnetic, is used to investigate surface mag- 
netism and to visualise domain walls. Although recent 
studies have attempted to measure energy dissipation be- 
tween an oscillating tip and a magnetic sample [9, 10], the 
dependence of the friction force on the tip's sliding velocity 
has not been considered yet apart from a work by C. Fusco 
et al. [8] which is extended by the present work to temper- 
atures r^O. The relative motion of the tip with respect 
to the surface can lead to the creation of spin waves which 



propagate inside the sample and dissipate energy, giving 
rise to magnetic friction. 

We will first present a simulation model and define mag- 
netic friction. The model contains classical Heisenberg 
spins located on a rigid lattice which interact by exchange 
interaction with each other. Analogous to the reading 
head of a hard disc or a MFM tip, an external fixed 
magnetic moment is moved across the substrate. Using 
Langevin dynamics and damping, it is possible to simu- 
late systems at finite temperatures. The main new results 
concern the temperature dependence of magnetic friction. 

Simulation model and friction definition. To 

simulate a solid magnetic monolayer (on a nonmagnetic 
substrate), we consider a two-dimensional rigid x Ly 
lattice of classical normalised dipole moments ("spins") 
Si = /x,j/^s, where fig denotes the material-dependent 
magnetic saturation moment (typically a few Bohr magne- 
tons). These spins, located at z = and with lattice spac- 
ing a, represent the magnetic moments of single atoms. 
They can change their orientation but not their absolute 
value, so that there are two degrees of freedom per spin. 
We use open boundary conditions. A constant point dipole 
Stip pointing in the z-direction and located at z = 2a rep- 
resents the magnetic tip. It is moved parallel to the surface 
with constant velocity v. 

This model has only magnetic degrees of freedom and 
thus focusscs on their contributions to friction. For a 
real tip one could expect that magnetic, just like nonmag- 
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netic [11,12] interactions might also lead to atomic stick- 
slip behaviour, and hence to phononic dissipation with 
a velocity-independent friction contribution as described 
by the Prandtl-Tomlinson model [13,14]. However, this 
requires a periodic potential between tip and substrate, 
that is strong enough compared to the elastic deformation 
energy to allow for multiple local potential energy minima. 
The magnetic tip-substrate interactions are unlikely to be 
strong enough. 

The Hamiltionan consists of two parts: 



H, — Tisub + Hsub-tip- 



(1) 



The first one represents the internal ferromagnetic short- 
range interaction within the substrate. The second one 
describes the long-range coupling between the substrate 
and the tip. 

The interaction between the substrate moments is mod- 
eled by the anisotropic classical Heisenberg model. 
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J > describes the ferromagnetic exchange interaction 
between two nearest neighbours, expressed by the an- 
gular brackets dz < quantifies the anisotropy. 
which prefers in-plane orientations of the spins. The 
dipole-dipole-interaction between the substrate spins is 
neglected, because it is much weaker than the exchange 
interaction. A quantitative comparison of our simulation 
results with the ones obtained in [8], where the dipole- 
dipole-interaction inside the substrate was taken into ac- 
count, justifies this approximation, which reduces simula- 
tion time enormously. 

The long-range interaction between substrate and tip is 
described by a dipole-dipole interaction term 



E3 (Si • e.i)(Stip • e^) — Si ■ Stip 


1=1 



(3) 



where Ri ~ |Ri| denotes the norm of the distance vector 
Ri = Ti — rtip, and ei its unit vector Ci = Ri/i?i. ri and 
rtip denote the position vectors of the substrate spins and 
the tip respectively, w quantifies the dipole-dipole cou- 
pling of the substrate and the tip. Note, however, that 
the results of the present study only depend on the com- 
bination w I Stip I, which is the true control parameter for 
the tip-substrate coupling. 

The proper equation of motion of the magnetic moments 
is the Landau-Lijshitz- Gilbert (LLG, [15]) equation, 



The first term on the right-hand side of eq. (4) describes 
the dissipationless precession of each spin in the effective 
field hi (to be specified below). The second term de- 
scribes the relaxation of the spin towards the direction 
of hi. 7 denotes the gyromagnetic ratio (for free electrons 
7 = 1.76086 X lO^^s^^T^^), and a is a phenomenological, 
dimcnsionless damping parameter. 

The effective field contains contributions from the tip 
and from the exchange interaction, as well as a thermally 
fluctuating term C,^ [17,18], 



(6) 



The stochastic, local and time-dependent vector Ci(0 '^X" 
presses a "Brownian rotation" , which is caused by the 
heat-bath connected to each magnetic moment. In our 
simulations this vector is realised by uncorrelated random 
numbers with a Gaussian distribution, which satisfy the 
relations 
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kBT5,,j5^,x5{t - t'), (8) 



where T is the temperature, (5ij- expresses that the noise 
at different lattice sites is uncorrelated, and (Sk,a refers to 
the absence of correlations among different coordinates. 

To find a quantity which expresses the friction occurring 
in the system, it is helpful to discuss energy transfers be- 
tween tip, substrate and heat-bath first. It is straightfor- 
ward to separate the time derivative of the system energy, 
eq. (1), into an explicit and an implicit one. 
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The explicit time dependence is exclusively due to the tip 
motion. The energy transfer between the tip and the sub- 
strate is expressed by the first term of eq. (9), which jus- 
tifies to call it the "pumping power" T'pump: 



on dH: 
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[S, X h, -I- a Si X (S, X h,)] , (4) At any instance, the substrate exerts a force 



'9'Waub-tip 



which is equivalent to the Gilbert equation of motion [16]: 



-S, X 



7 dt 



(5) 



on the tip. Due to Newton's third law, T^pump is the work 
per unit time done by the tip on the substrate. Its time 
and thermal average ('Ppump) is the average rate at which 
energy is pumped into the spin system. In a steady state 
it must be equal to the average dissipation rate, i.e. to 
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Figure 1: A single spin in a magnetic field rotating with angular 
velocity <jJ, is dragged along with a phase shift Sip and 
an out of plane component SO. 
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the net energy transferred to the heat bath per unit time 
due to spin relaxation. The magnetic friction force can 
therefore be calculated by 



F : 



pump / 



(11) 



The second term of eq. (9) describes the energy transfer 
between the spin system and the heat bath. Inserting eq. 



(5) into Pdiss = - E 
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(12) 



The first term, "Pthorm containing describes, how much 
energy is transferred into the spin system due to the ther- 
mal perturbation by the heat bath. The second term, 
7'rciax proportional to the damping constant a, describes 
the rate of energy transfer into the heat bath due to the 
relaxation of the spins. 

At T = 0, Vtherm IS zero. The spins are only perturbed 
by the external pumping at u 0. Then 
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(T = 0), 
(13) 



where for the last transformation we used eq. (4) in order 
to show explicitly the relationship between dissipation rate 
and misalignment between spins and local fields. 

It will be instructive to compare the magnetic substrate 
scanned by a dipolar tip with a much simpler system, in 



which the substrate is replaced by a single spin S sub- 
jected to an external field h{t) that rotates in the plane 
perpendicular to a constant angular velocity a; (replac- 
ing the tip velocity). It is straight forward to obtain the 
steady state solution for T = 0, where in the co- rotating 
frame S is at rest. S lags behind h/h by an angle 6tp and 
gets a component Sd in u;-direction (cf. fig. 1), which are 
in first order given by 
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Inserting this into the (iV=l)-case of eq. (13) yields a dis- 



sipation rate of T^diss 



which corresponds to a 



"viscous" friction F ~ Vdiss/^ oc auj. 

It is instructive to give a simple physical explanation 
for eq. (14). instead of presenting the general solution, 
which can be found in [19]. Two timescales exist in the 
system, which can be readily obtained from eq. (4); first, 
the inverse Lamor frequency TLamor = (1 + Q!^)Ms/7^i ^-nd 
second, the relaxation time Treiax=TLamor/Q;. They govern 
the time evolution of 5ip and S0. In leading order. 
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(15) 
(16) 



The first equation describes, how S9 would increase by pre- 
cession of the spin around the direction of the field, which 
is counteracted by relaxation back towards the equator. 
The second equation describes that without relaxation 
into the field direction, Sp would increase with velocity 
uj minus the azimuthal component of the precession veloc- 
ity, which is in leading order proportional to S9. Setting 
the left hand sides to zero in the steady state, immediately 
gives the solution (14). 

In the (T = 0)-study [8], the time average {Vrciax)/v was 
used to calculate the friction force. As pointed out above, 
this quantity agrees with (11) in the steady state. For 
finite temperatures, however, (11) is numerically better 
behaved than {Pieiax)/v- The reason is the following: 

For T 7^ 0, the spins are also thermally agitated, even 
without external pumping, when the dissipation rate 'Pdiss 
vanishes. This shows that the two terms Ptherm and Preiax 
largely compensate each other, and only their difference is 
the dissipation rate we are interested in. This fact makes 
it difficult to evaluate (12) and is the reason why we prefer 
to work with (11) as the definition of the friction force. 

We have also analyzed the fiuctuations of the friction 
force (11). The power spectrum has a distinct peak at 
frequency v/a, which means that the dominant temporal 
fiuctuations are due to the lattice periodicity with lattice 
constant a. A more complete investigation of the fluctu- 
ations, which should also take into account, how thermal 
positional fiuctuations influence the friction force, remains 
to be done. 
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Technical remarks. — Because of the vector pro- 
duct in eq. (4), the noise Ci enters in a multipUcative 
way, caUing for special attention to the interpretation of 
this stochastic differential equation (Stratonovich vs. Ito 
sense). The physical origin of the noise renders it gcneri- 
cally coloured and thus selects the Stratonovich interpre- 
tation as the appropriate one (Wong-Zakai theorem [20]), 
in which its appearance as white noise is an idealisation. 
Accordingly we employ the Hcun integration scheme [21]. 
After each time step the spins are rescaled so that their 
length remains unchanged. 

To get meaningful results, it is of prime importance to 
reach a steady state. The initial configuration turned out 
to be a crucial factor for achieving this within acceptable 
computing time. Therefore a long initialisation run is per- 
formed, before the tip motion starts. 

Moreover, the system size is another limiting factor. In 
order to avoid that the tip reaches the system boundary 
before the steady state is reached, we use a "conveyor belt 
method" allowing to do the simulation in the comoving 
frame of the tip. The tip is placed in a central point above 
the substrate plane, e.g. at {{Lx + l)/2,{Ly + l)/2). After 
an equilibration time, the tip starts to move with fixed 
velocity in x-dircction. When it passed exactly one lattice 
constant a, the front line (at x = L^) is duplicated and 
added to the lattice at x = Lj; + I. Simultaneously, the 
line at the opposite boundary of the system (at x ~ 1) is 
deleted, so that the simulation cell is of fixed size and con- 
sists of the Lx X Ly spins centered around the tip position, 
with open boundaries. Note that this is different from pe- 
riodic boundary conditions, because the spin configuration 
deleted at one side is different from the one added at the 
opposite side. We compared the results obtained for a 
small system in the co-moving frame of the tip with some 
runs for a system that was long enough in x-direction that 
the steady state could be reached in the rest frame of the 
sample, and confirmed that the same friction results and 
steady state properties could be obtained with drastically 
reduced computation time. 

It is convenient to rewrite the equations of motion in 
natural units. An energy unit is prescribed by the ex- 
change energy J of two magnetic moments. It rescales the 
energy related parameters dz and w as well as the simu- 
lated temperature. 



knT' 



J 



(17) 



The rescaled time further depends on the material con- 
stants and 7, 



t' ^ — t. 

fJ.s 



(18) 



A length scale is given by the lattice constant a, so a nat- 
ural velocity for the system can be defined. 




0.3, v = 0.05 
0.3, v = 0.10 
0.5, v = 0.10 
tanh(0.75 x) 



Figure 2: Local magnetisation (x-component) at T—Q along 
the lattice axes in a^-direction which are closest to the tip (i.e. 
at 2/=±0.5). Analogous to a domain wall one finds a tanhs- 
profile. Depending on the damping constant a and the velocity 
V, its zero is shifted backwards from the tip position by a value 
xo~—0.88av, as shown in the inset. 

From the natural length and energy, the natural force re- 
sults to: 



F' 



J 



-F 



(20) 



From now on, all variables are given in natural units, and 
we will drop the primes for simplicity. The typical exten- 
sion of the simulated lattices is 50 x 30, which was checked 
to be big enough to exclude finite size-effects. For the tip 
coupling, we chose large values (e.g. wStip = (0,0,-10)), 
to get a large effective field on the substrate. Usually it 
is assumed that the dipole-dipole coupling constant has a 
value of about w = 0.01, which means that the magnetic 
moment of the tip is chosen a factor of 1000 times larger 
than the individual substrate moments. The anisotropy 
constant is set to dz — —0.1 in all simulations. The damp- 
ing constant a is varied from 0.1 to the quite large value 1. 
At finite temperatures typically 50 simulation runs with 
different random number seeds are performed to get reli- 
able ensemble averages. 

Simulation results. — In [8] it was found that the 
magnetic friction force depends linearly on the scanning 
velocity V and the damping constant a for small velocities 
(u<0.3). For higher velocities the friction force reaches a 
maximum and then decreases. In this work we focus on 
the low-velocity regime with the intention to shed more 
light on the friction mechanism and its temperature de- 
pendence.^ 

Adiabatic approximation at T— . If we assume the 
field h for each spin to vary slowly enough to allow the so- 
lution (14) to be attained as adiabatic approximation, the 
linear dependence F (x av from [8] follows immediately: 
At every point, the temporal change of the direction of 
h, defining a local to for (14), is proportional to v. We 



, Ms 
V ~ V. 



(19) 



^It should be noted that the smallest tip velocities we simulated, 
are of the order of 10~^ (aj'y / ^s), which is still fast compared to 
typical velocities in friction force microscopy experiments. 
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Figure 3: Left: Magnetisation profiles as in fig. 2 for several temperatures with wS'tip = 10, a = 0.5 and v = 0.01. Middle and 
right: Absolute value of the average magnetisation as a function of the distance r from {x,y) = (0,0), directly underneath the 
tip. For small temperatures (upper two curves) it decreases with a power law (cf. double-logarithmic plot, middle), for high 
temperatures (lower two curves) it decreases exponentially (cf. semi-logarithmic plot, right). 



confirmed the validity of the adiabatic approximation nu- 
merically by decomposing S ~ h/h with respect to the 
local basis vectors dt{h/h), h, and their cross-product, all 
of them appropriately normalized. In other words, we ex- 
tracted 5ip and SO directly and found them in excehent 
agreement with (14). 

The lag of S with respect to h manifests itself also 
macroscopically in the magnetisation field as we will show 
now. The tip-dipole is strong enough to align the substrate 
spins to nearly cylindrical symmetry. Since the anisotropy 
is chosen to generate an easy plane (c?z<0), spins far away 
from the tip try to lie in the xy-plane, while close to the 
tip they tilt into the z-direction. This is displayed in 
fig. 2 where the z-component of the local magnetisation is 
shown along a line in x-direction for a fixed y-coordinate. 
Remarkably, these magnetisation profiles for different val- 
ues of V and a collapse onto a unique curve, if they are 
shifted by corresponding offsets xq with respect to the tip 
position. As expected from (14), the magnetisation profile 
stays behind the tip by a (^/-dependent) shift xq oc av (cf. 
inset of fig. 2). 

Friction at T>0. With increasing temperature the 
magnetisation induced by the tip becomes smaller, as 
shown in fig. 3. One can distinguish a low tempera- 
ture regime, where the local magnetisation decreases alge- 
braically with the distance from the tip, and a high tem- 
perature regime, where it decreases exponentially. The 
transition between these regimes happens around T w 0.7. 

For all temperatures the friction force F turns out to be 
proportional to the velocity (up to v ~ 0.3), as for T=0, 
with a temperature dependent friction coefficient F/v. 

The two temperature regimes manifest themselves also 
here, as shown in fig. 4: For low temperatures the fric- 
tion coefficients depend nearly linearly on a, reflecting the 
T=0 behaviour. Towards the high temperature regime, 
however, the a-dependence vanishes, and all friction coef- 
ficients merge into a single exponentially decreasing func- 
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wS^i = 10, a=0.5 
wi'up = 10, a=0.7 
vvSti = 5, a=0.5 
wStip= 5, a=0.7 




0.2 0.4 0.6 



0.8 



1.2 



1.4 1.6 



Figure 4: Friction coefficients for different a, wStip and ksT. 
One can distinguish between a low temperature regime, where 
the friction coefficient depends on a but not on wStip, and a 
high temperature regime, where it depends on wSup but not 
on a. 



tion. 

The low temperature behaviour can be understood es- 
sentially along the lines worked out for T=0, as result of 
a delayed, deterministic response (precession and relax- 
ation) to the time dependent tip field. At high tempera- 
tures, however, friction results from the ability of the tip to 
propagate partial order through the thermally disorderd 
medium. The magnetisation pattern in the wake of the tip 
no longer adapts adiabatically to the dwindling influence 
of the tip, but decays due to thermal disorder. Then, a ris- 
ing temperature lets the ordered area around the tip shrink 
which leads to the exponential decrease of the friction coef- 
ficient. However, it increases with the tip strength, wStip, 
as stronger order can be temporarily forced upon the re- 
gion around the tip. By contrast, the tip strength looses 
its influence on friction in the limit T — > 0, because the 
substrate spins are maximally polarised in the tip field . 

This picture of the two temperature regimes is sup- 
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Figure 5; Distance Xq by which the magnetisation pattern lags 
behind the tip is proportional to v for all temperatures. The 
proportionality constant depends on a only in the low temper- 
ature regime. 

ported by the distance xq, by which the magnetisation 
pattern lags behind the tip. It is proportional to v for all 
temperatures, but a-dependent only in the low tempera- 
ture regime, cf. fig. 5. 

Conclusion and outlook. In this work, we could 
explain the low-velocity, zero-temperature findings from 
[8] , namely that the magnetic friction force in the Heisen- 
berg model has a linear velocity dependence with a coeffi- 
cient proportional to the damping constant a. In the spin 
polarisation cloud dragged along with the tip, each sub- 
strate spin follows the local field with a lag proportional to 
the frequency of the field change and to a. Moreover, the 
magnetisation pattern around the tip gets distorted due to 
precession. These effects directly give rise to the observed 
magnetic friction and could be evaluated quantitatively by 
means of a single spin model. 

Second, for the first time the temperature dependence 
of magnetic friction in the Heisenberg model was investi- 
gated in the framework of Landau-Lifshitz-Gilbert (LLG) 
dynamics with a stochastic contribution to the magnetic 
field. Two regimes were found, which can be charac- 
terised by their different relaxation behaviour. While in 
the low-temperature regime the response of the system on 
the perturbation due to the moving tip is dominated by 
the deterministic precession and relaxation terms in the 
LLG equation, thermal perturbations competing with the 
one caused by the moving tip are essential in the high- 
temperature regime. This explains, why magnetic friction 
depends on a but not on wSfip for low temperatures, while 
it depends on wSup but not on a for high temperatures 
where it decreases exponentially with T. 

Important extensions of the present investigation in- 
clude the effects of a tip magnetisation pointing in a differ- 
ent than the z-direction, of the strength and sign of spin 
anisotropy, dz, or of the thickness of the magnetic layer. 
Both, spin anisotropy and lattice dimension will be crucial 
for the critical behaviour, as well as for the critical tem- 



perature itself. Studies dealing with these quantities are 
already in progress and will be reported in a future work. 
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